Coulomb correlations and coherent charge tunneling in mesoscopic coupled rings 



C. M. Canali 1 , W. Stephan 2 , L. Y. Gorelik 1 , R. I. Shekhter 1 and M. Jonson 1 
1 Dept. of Applied Physics, Chalmers University of Technology and Goteborg University, S-412 96 Goteborg, Sweden 
2 Max- Planck- Institut fur Physik komplexer Systeme, Bayreuther Str. 40, D-01187 Dresden, Germany 

(1 October 1997) 

We study the effect of a strong electron-electron (e-e) interaction in a system of two concentric 
one- dimensional rings with incommensurate areas A\ and A2, coupled by a tunnel amplitude. For 
noninteracting particles the magnetic moment (persistent current) m of the many-body ground state 
and first excited states is an irregular function of the external magnetic field. In contrast, we show 
that when strong e-e interactions are present the magnetic field dependence of m becomes periodic. 
In such a strongly correlated system disorder can only be caused by inter-ring charge fluctuations, 
controllable by a gate voltage. The oscillation period of m is tx 1/(j4i + A%) if fluctuations are 
suppressed. Coherent inter-ring tunneling doubles the period when charge fluctuations are allowed. 

PACS numbers: 71.27.+a, 73.20.Dx, 72.15.Nj 



Persistent currents in quasi onedimensional mesoscopic 
rings jjj are a phase sensitive probe of the conductor 
equilibrium properties. The intriguing role played by e-e 
interactions in such a phase coherent phenomenon con- 
tinues to be the source of intensive investigation. On 
one hand theoretical work || has shown that interaction 
does not influence the persistent currents in clean single- 
channel rings, in agreement with experiments on a single 
ballistic semiconductor ring with few channels On 
the other hand many other theoretical studies 0] have 
explored the possibility of linking the e-e interaction to 
the puzzlingly large current amplitudes found experimen- 
tally in disordered metallic (i.e. many-channel) rings 
Such amplitudes cannot be accounted for within a single- 
particle theory. No general consensus has been found yet 
on the effects of interaction in these systems. 

As an example where e-e interaction is expected to 
play a distinctive role and cause interesting effects, we 
consider here a system of two concentric ballistic rings, 
coupled by a tunnel amplitude. The areas enclosed by 
the two rings are A\ and A2 respectively, and a per- 
pendicular magnetic field B is present. Our aim is to 
investigate how strong interactions affect the many-body 
ground-state and the first excited states as a function of 
the applied magnetic field. We will show that the interac- 
tion creates a spatially ordered rigid structure, where the 
two rings are tied together and charge fluctuations tend 
to be suppressed. Consequently the flux dependence of 
the lowest energy levels becomes periodic with period 
cx \/{Ai + A2), in contrast to the irregular aperiodic 
behavior of the free-particle case. By means of a gate 
voltage one can induce the coherent transfer of a single 
electronic charge from one ring to the other. As a result 
of the coherent charge fluctuations between the rings, the 
oscillation period of the lowest levels is doubled. Charge 
transfer induced effects in different mesoscopic samples 
have been recently studied by Buttiker and Stafford || . 

We choose a lattice model of two concentric one- 
dimensional rings with AT sites per ring (see Fig. 1), de- 
scribed by the following second quantized Hamiltonian 
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The operators cj and c i create and destroy a spin- 
less electron at site i of ring p. (i, j — 1, . . . , Af are site 
indexes and p,p' — 1,2 are ring indexes). The term pro- 
portional to t describes the intra-ring hopping, affected 
by the magnetic field through the phases ip p — BA p /4>o, 
where cf>o — ch/e is the flux quantum. The site-dependent 
amplitude £12(2) represents inter-ring tunneling. The 
last term is the translationally invariant e-e interaction. 
The number of electrons in each ring, N p , is not a good 
quantum number if the rings are coupled. The expec- 
tation values (N p ) will however satisfy the constraint 
(Ni) + (N2) = N; N is the total number of electrons 
in the system. 

We first discuss the case of two decoupled [<i2(i) = 0] 
noninteracting rings. The single-particle spectrum of 
each ring is given by e p (fcj) = — 2tcos[(ki + 2irtp p /Af)] = 
-2icos[(fc; + 2nH a p /Af)], with fe, = i = 

0, ±1, . . . ± Af/2. Here we have introduced two indepen- 
dent, dimensionless parameters, H (magnetic field) and 
o p (area of ring p) . When H — 1 and a p = 1 the flux en- 
closed by ring p is <fio . As H varies the levels move up and 
down; at some values of H they can cross. As a result, 
the particle at the Fermi energy can switch to an empty 
state of different "velocity", causing a cusp in the mag- 
netic field dependence of the ground state energy, Eq (H) . 
Cusps in Eq(H) will cause jumps in the total magnetic 
moment of the ground state, m(H) = dE$(H) / dH . 
When N is even one can see that, in the ground-state, 
each ring of the system contains always N/2 particles. 
When flux-induced level crossing occurs, the particle at 
the Fermi energy will change place to an empty state of 
the same ring with opposite velocity. Therefore m(H) 
will consist of the superposition of the contributions of 
two independent rings, each one oscillating with its own 
period l/a p as a function of H . If the areas of the two 
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rings are incommensurate, - throughout this work we 
have considered the values a\ = [\/(5) — l]/2 w 0.618 
and a2 = 1.5 - the resulting total m(H) is a bound 
but aperiodic function of H, with sharp jumps occur- 
ring at multiple of l/a p , p — 1,2. On the other hand 
when AT is odd the relevant level crossings involve al- 
ways two states belonging to different rings: an occupied 
state of the more populated ring is depleted in favor of 
an empty state of the other ring. In this case the ape- 
riodic m(H) contains jumps at values of H which are 
multiples of l/(oi ±02), corresponding to level crossings 
between states of different rings moving in opposite or 
equal direction respectively. 

To first order in perturbation theory a small tunnel- 
ing amplitude |tia(i)l << A = 8i/A/ (here A is the 
mean level spacing for each single ring) will couple states 
of different rings. If ii2(t) is also site independent, the 
spectrum will hardly be modified and, in particular, the 
jumps in m{H) will remain sharp, since they always in- 
volve crossing between levels of different momenta which 
cannot be coupled by a homogeneous hopping. However, 
if t\2(i) is not homogeneous, states of different momenta 
can be coupled, causing a rounding of the sharp jumps in 
m(H). If I tx2 I ~ A the energy levels of the two rings 
will be strongly hybridized and m(H) will depend on the 
specific form of fi2(*)- 

In order to see how strong e-e interactions modify this 
scenario we have resorted to exact numerical calculations 
on small systems, using the Lanczos algorithm to com- 
pute the many-body ground state and the first low-lying 
excitations. We have considered systems with Af = 11 
sites per ring and N = 4, 5, 6, 7, 8 particles. We will first 
discuss the case of a short range interaction, coupling 
two nearest-neighbor sites (i,p) and (i + \,p) in the same 
ring with matrix element V\\, and opposite sites (i, 1) and 
(i, 2) in different rings with matrix element V\2- The rel- 
evant energy parameter, indicating the relative strength 
of the interaction, is V pp >/t. Plots of the magnetic mo- 
ment for systems of N = 6, 7, 8 particles interacting with 
Vn = 6.5 1, V12 — 6.3 1 are shown in Fig. 2. The ef- 
fect of the interaction is visibly spectacular, particularly 
when iV is even: the aperiodic -ff-dependence of m(H) 
of the noninteracting case is replaced by a perfectly peri- 
odic function, oscillating with period Hq — l/{a% +02), 
exactly as we would have for a single ring with area 
fltot = (fli + 0,2)- The periodic m(H) develops smoothly 
upon increasing the interaction and once it is reached it 
becomes independent of the interaction strength. The 
amplitude of m(H) in the strongly interacting regime is 
approximately one half of the amplitude of the nonin- 
teracting case. When N is even, the persistent current 
1(H) = —cm(H) is diamagnetic, opposite to the single- 
ring case. As shown in Fig. 2(c), an inter-ring hopping 
t\2(i) of any kind does not spoil this result, even when 
its magnitude is comparable to the intra-ring hopping 
t. Morover, not only is the ground-state energy a peri- 
odic function of H, but so are the first low-lying excited 
states. In particular, the low-energy sector of the many- 



body spectrum is qualitatively identical to the single- 
particle spectrum in one ring. We can also evaluate the 
particle number in each ring for each many-body state 
by computing the corresponding wavefunction. The cal- 
culation shows that in the strongly interacting regime all 
periodic eigenstates have Ni = N2 = N/2 for any value 
of H when the rings are decoupled. Inter-ring charge re- 
distributions, causing _ff-dependent deviations from the 
(A^x = N/2,N2 — N/2) configuration, are suppressed 
here. Charge fluctuations occur however in the higher 
aperiodic states; these are separated from the ordered 
low-energy sector by an energy gap that scales with 
Vppi ft. A non-zero £12(2) ~ A << V pp i cannot connect 
these two sectors of the spectrum. Therefore the low- 
lying states remain periodic with {N\) » (A^) ~ N/2, 
even when the rings are coupled. 

Strong intra-ring interaction is known to generate 
Wigner-crystal-like ground states in a single ring [Q. It 
has also been shown that inter-ring interaction causes a 
nondissipative Coulomb drag H that affects qualitatively 
the persistent currents in each ring || . Our exact results 
suggest that strong inter-ring interaction glues together 
the two stiff Wigner crystals created by the intra-ring in- 
teraction, resulting in a rigid structure that rotates as a 
single solid body under the effect of the magnetic field. 
Based on this picture, we can easily see why the oscil- 
lation period of the total magnetic moment of such a 
structure is l/(di +02)- Indeed the Hamiltonian of the 
two frozen crystals glued together, each containing N/2 
particles, is found to be 
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where —ih-^g is the canonical momentum conjugated to 
6, the rotation angle of the solid around an axis pass- 
ing through the center of the rings; X — m e ^r(^4i + A 2 ) 
is the moment of inertia. In zero flux, the eigenstates 
^n(O) = exp(in9) must satisfy the boundary condition 
$r„(0) = $ n (6 + 2w/(N/2)) due to the ordered crys- 
tal structure. This immediately implies the announced 
l/(ai + CI2) periodicity as a function of H. 

The m(H) for odd N [see Fig 2(b)] becomes also pe- 
riodic in presence of a strong interaction, with the same 
period Hq. One significant difference with respect to even 
A^ is the shift of m(H) by Ho/ 2: the system is now para- 
magnetic. 

We have also studied a more realistic unscreened 
long range interaction coupling all the sites with ma- 
trix element V^(|rj iP — Tj >p r\) = V c Ri/\r itP — r,- 

V cy /a\/ a p + a p i — 2^/cip a p ' cos(27r(j — j)/Af) 
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dimensionless interaction constant is a = V c R\/Tiv-p = 
(V^/t)[A/"/47rsin(A;p)]. In the strong interaction limit 
(a > 5) the Hamiltonian Eq. [1] has an unphysical ground 
state where all the particles are located in the outer ring. 
Obviously a realistic Hamiltonian should contain the con- 
tribution of a positive background preventing the deple- 
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tion of the inner ring. This is usually done by replac- 
ing 7ii )P by rii.p — K in Eq. [1], where K represents the 
average positive jellium which equals the average elec- 
tron density. Such a replacement is formally equiva- 
lent to add in the Hamiltonian Eq. [1] a potential dif- 
ference between the rings, H g = J2i P e p n i,p- Since we 
also want to consider the effect of a real external gate 
voltage V g — ei — £2 we will assume that the background 
is incorporated in V g , which we will take as a parame- 
ter. Typically a voltage difference V g » — V c /2 ensures 
an almost equal population in the two rings. In Fig. 3 
we plot m(H) for N = 5, 6 particles interacting with a 
Coulomb force of intensity V c = 10 1, for two values of 
V g . Let us consider first N — 6 and t\2 = (black line 
in Fig. 3(a)(b)). For V g = -0.75 V c , m(H) is perfectly 
periodic with the expected period 1 / (ai + a 2 ). On the 
other hand imperfections in the periodic pattern appear 
for V g = —0.5 V c . Note in particular the small glitches 
at H = 0.3, 1 and the splitting of the peak at H = 3 in 
Fig. 3(a). In fact the calculation of the ring occupancy 
shows that while for Vg = — 0.75 V c the electronic configu- 
ration (JVj. = 3, N 2 = 3) is stable for all values of H, when 
V g = 0.5 V c the two configurations (Ni = 3, N 2 = 3) and 
(Ni = 2,N 2 = 4) alternate in the ground state. The 
aperiodic jumps in Fig. 3(a) correspond to a switch from 
one configuration to another as H varies. It is instruc- 
tive to look at the H dependence of the lowest four energy 
levels, plotted in Fig. 4, and their corresponding ring oc- 
cupancy. Fig. 4(a) shows that when V g = —0.5 V c these 
four levels are often two by two degenerate, the degen- 
eracy being caused by the presence of the two electronic 
configurations (Ni = 3, N 2 = 3) and (Ni = 2,N 2 = 4) 
mentioned above. V g is not yet strong enough to en- 
force a fixed electronic configuration in the two rings. 
As a result the levels are not completely periodic and 
they are sensitive to inter-ring tunneling- see grey line 
of Fig 3(a). However if we increase V g up to 0.75 V c , 
the configuration (Ni = 3, N 2 = 3) is locked in the first 
three states. Charge redistribution becomes possible only 
starting from the fourth level, with the appearance of 
the (N± = 2,N 2 = 4) configuration at some values of H. 
The important consequence is the ordering of the first 
three levels, which is robust against tunneling between 
the rings, as shown in Fig. 3(b) and Fig. 4(c). 

The example of 6 particles is representative of what 
should happen in a real mesoscopic system with N% ~ 
N 2 ~ 50, where gate- and flux-induced charge redistri- 
bution between the rings amounts to small changes in 
their electron concentration. If the gate voltage enforces 
an almost equal population and prevents charge fluctua- 
tions between the rings as H varies, then the e-e Coulomb 
interaction can fully develop the ordering effect. 

The presence, at some values of Vg, of two quasi- 
degenerate states that differ by one single electron being 
moved from one ring to the other gives rise to the possi- 
bility of coherent tunneling Q in presence of strong cor- 
relation. This effect is best illustrated for odd N when 
the swapping of one electron between the rings gener- 



ates the two configurations (N% = N 2 — ^^-)- 

In this case a homogeneous ti 2 coupling the two rings 
allows the swapping electron to occupy an hybridized or- 
bital involving two opposite sites of the two rings with 
equal probability. Since the particles are indistinguish- 
able all of them can be coherently involved in this pro- 
cess. This means that all the N particles will spend half 
of their time on either ring, that is (N±) w (N 2 ) w N/2, 
which is equivalent to having N particles on a single ring 
of area (a\ + a 2 )/2. The period of the persistent cur- 
rent of such "dynamically-glued" Wigner crystals will be 
2H = 2/(oi + a 2 ), twice as large as the previous case, 
where charge fluctuations were totally suppressed. The 
numerical results for TV = 5 shown in Fig. 3(c) confirms 
fully this prediction: the presence of a homogeneous t\ 2 
coupling coherently the quasi-degenerate configurations 
(Ni = 2,N 2 = 3) and (Ni = 2, N 2 = 3) doubles the oscil- 
lation period of m(H). Furthermore also the amplitude 
of m(H) is factor of 2 larger than when i 12 = 0. A sim- 
ilar result holds approximately for N — 6, see Fig. 3(b), 
but now the oscillating period must be generalized to 
N/((Ni)ai + (N 2 )a 2 ). For N — 6 the computed averaged 
occupations are (Ni) w 2.5, (N 2 ) w 3.5 and the resulting 
period is a little smaller than 2/(ai +02). When N >> 1 
we have {Ni)/N w (N 2 )/N w 1/2 and we obtain the re- 
sult that holds exactly for N = 5. 

Our calculations show that spectrum ordering is al- 
ready well developed for a ss 5. Experimentally, GaAs- 
Al/GaAs heterostructures with a > 1 are available. Thus 
the effects that we have found might be seen in a sys- 
tem of two semiconductor rings. Our results might also 
be connected to the halving of the oscillation period of 
the Aharonov-Bohm effect in Quantum Hall edge states 
around a semiconductor antidot Jiof . The experiment in 
Ref. |1(J showed that when just two edge states encircled 
the antidot at approximately the same radius, two sets of 
AB oscillations seemed to lock together, but out of phase, 
giving an apparent (h/e)/2 oscillation. This seems to 
correspond to a double-ring system with A\ — A 2 , with 
periodicity proportional to l/2Ai. 

In conclusion, we have shown that e-e interaction or- 
ders the flux dependence of the low-energy spectrum of 
a two-ring system. When strong e-e correlations are 
present, disorder can only be induced by inter-ring charge 
redistribution controllable through an external gate. 

We thank I. V. Krive for useful discussions. Financial 
support from NFR and TFR is gratefully acknowledged. 



[1] M. Biittiker, Y. Imry, and R. Landauer, Phys. Lett. 96A, 
365 (1983); M. Biittiker, Phys. Rev. B 32, 1846 (1985). 

[2] A. J. Legget, in Granular Nanoelectronics, edited by 
D. K. Ferry (plenum Press, New York, 1991), p. 297. 
D. Loss, Phys. Rev. Lett. 69, 343 (1992). 



3 



[3] D. Mailly, C. Chapelier, and A. Benoit, Phys. Rev. Lett. 
70, 2024 (1993). 

[4] V. Ambegaokar and U. Eckern, Phys. Rev. Lett. 65, 381 
(1990); A. Schmid, ibid. 66, 80 (1991); B. L. Altshuler, 
Y .Gefen, and Y. Imry, ibid. 66, 88 (1991); M. Abraham 
and R. Berkovits, ibid. 70, 1509 (1993). 

[5] L. P. Levy, G. Dolan, J. Dunsmuir, and H. Bouchiat, 
Phys. Rev. Lett. 64, 2074 (1990); V. Chandrasekhar, 
R. A. Webb, M. J. Brady, M. B. Ketchen, W. J. Galager 
and A. Kleinsasser, ibid. 67, 3578 (1991). 

[6] M. Biittiker and C. A. Stafford, Phys. Rev. Lett. 76, 495 
(1996). 

[7] H. J. Schultz, Phys. Rev. Lett. 71, 1864 (1993); 

I. V. Krive et al, Phys. Rev. B 52 16451 (1995). 
[8] A. G. Rojo and G. D. Mahan, Phys. Rev. Lett. 68, 2074 

(1992). 

[9] T. V. Shahbazyan and S. E. Ulloa, unpublished. 
[10] C. J. B. Ford et al. Phys. Rev. B 49, 17456 (1994). 



FIG. 2. Magnetic moment vs H for the sample of Fig. 1 
with Af = 11 sites per ring. The normalized areas are 
«i = [\ft5) - i]/ 2 « 0.618 and a 2 = 1.5 respectively. The 
total number of particles is: (a) N = 6, (b) N = 7 
and (c) N — 8. Aperiodic dashed curves: decoupled rings 
[ti 2 (i) = 0], noninteracting particles. Periodic black lines: 
ti2{i) = 0, particles interacting via a short range force; the 
oscillation period, is Ho = l/(ai + a 2 ) ~ 0.47. Periodic gray 
line in (c): £12(1) = — t/2*cos 2 (3ni/2J\f), interacting particles. 
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FIG. 3. Magnetic moment vs H for particles interacting 
with a Coulomb potential of intensity V c = 10 t. V g is the 
voltage between the rings. Black lines: £12(2) = 0; Gray lines: 
ti 2 (i) = -t/2. (a) N = 6 particles, V g = -0.514. (b) N = 6, 
(c) N = 5, V g = -0.5 V c . 



FIG. 1. Schematic picture of two concentric rings with en- 
closed areas Ai and A 2 . The black dots are the lattice sites. 
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FIG. 4. Energy vs H of the ground state and the first three 
excited states for a system of N — 6 particles, interacting with 
a Coulomb potential of intensity V c = Wt. V g is the voltage 
between the rings, (a) £12 = 0, V g = 0.5V C . (b) ti 2 = 0, 
V g = -0.75 V c . (c) h 2 (i) = -t/2, V g = -0.75 V c . 
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